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Abstract 

Type I X-ray bursts from low-mass X-ray binaries result from a thermonuclear runaway in the 
material accreted onto the neutron star. Although typical recurrence times are a few hours, consistent 
with theoretical ignition model predictions, there are also observations of bursts occurring as promptly 
as ten minutes or less after the previous event. We present a comprehensive assessment of this 
phenomenon using a catalog of 3387 bursts observed with the BeppoSAX/WFCs and RXTE/PCA 
X-ray instruments. This catalog contains 136 bursts with recurrence times of less than one hour, that 
come in multiples of up to four events, from 15 sources. Short recurrence times are not observed from 
so-called ultra-compact binaries, indicating that hydrogen burning processes play a crucial role. As far 
as the neutron star spin frequency is known, these sources all spin fast at over 500 Hz; the rotationally 
induced mixing may explain burst recurrence times of the order of 10 min. Short recurrence time 
bursts generally occur at all mass accretion rates where normal bursts are observed, but for individual 
sources the short recurrence times may be restricted to a smaller interval of accretion rate. The fraction 
of such bursts is roughly 30%. We also report the shortest known recurrence time of 3.8 minutes. 
Subject headings: accretion, accretion disks — methods: observational — stars: neutron — X-rays: 
binaries — X-rays: bursts 



1. INTRODUCTION 

Type I X-ray bursts are thought to result from 
thermonuclear flashes of hydrogen and/or heli um in 
the envelope of neutron stars (iW oosley & Taam 119761 : 
iMaraschi fc Cavalierdfl977l : iLamb fc Lambl ll9781. This 
material is accreted through Roche-lobe overflow from 
a lower-mass companion star (low-mass X-ray binary, 
LMXB). Current one-dimensional models successfully 
explain burst features such as the peak fl ux, the fluence, 
decay time and recurre nce time (e.g., 'Woosley et al.' 
2001 IHegereFaT] [20071: see I Wallace fc Woo slev 1981; 



Fuiimoto et all 119811: iFushiki fc LambI 119871 for earlier 



work). During the flash, over 90% of the accreted hydro- 
gen and heliu m is expected to burn to carbon and heavier 
elements (e.g.. lWoosley et al]l2004n . For the next flash to 
occur, a fresh layer of hydrogen/helium must first be ac- 
creted. At typical accretion rates of up to approximately 
10~* M0yr~^ this takes at least a few hours. 

X-ray bursts have been observed sin ce the 1970's 
(jGrindlav et al.l 119761 : iBeHan et al.l Il976f ) from approx- 
imately 90 sources in our Galax y, with recurrence 
times of hours up to days ( e.g.. iLewin et all 119931 : 
iStrohmaver fc Bildstenl I2OO60 . ILewin et all (|1976D re- 
ported the detection with SAS-3 of three bursts that 
were separated by only 17 and 4 minutes. These bursts 
originated from a crowded region and, therefore, source 
confusion cannot be ruled out. In the 1980's similar re- 
currence times as short as 10 min utes were observed from 
both 4U 1608-522 with Hakucho (iMurakami e t al. 1980) 
and from EXO 0748-676 with EXOSAT (Ic'ottwald et al. 
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119861 Il987a[ ). This rare phenomenon implies that hy- 
drogen and helium is left over somewhere on the star 
after the initial flash, because the recurrence time is too 
short to accrete enough fuel for the subsequent burst (s). 
This is at odds with the current models, that predict an 
almost complete burning of the available hydrogen and 
hel ium on the e ntire s tar surface. 

Boi rin et all (|2007f) analyzed 158 hours of XMM- 
Newton observations of EXO 0748-676, which revealed 
short recurrence time bursts in groups of two (doubles) 
and three (triples). This relatively large burst sample 
revealed that on average bursts with a short recurrence 
time (8 to 20 minutes) are less bright and energetic than 
bursts with 'normal' recurrence times (over 2 hours). The 
fit of a black body model to the burst spectrum shows 
a lower peak temperature, while the emitting area is the 
same. The profiles of short recurrence time bursts seem- 
ingly lack the long 50 s to 100 s tail caused by rp-process 
burning, which indicates that the burst fuel contains less 
hydrogen. After a double or triple it takes on average 
more time before another burst occurs, suggesting a more 
c omplete burning of the available fuel. 

iGallowav et al.l ()2008[ ) showed that there are more 
sources that show this behavior, with bursts occurring 
in groups of up to four bursts, and with recurrence times 
as short 6.4 minutes. The short recurrence times were ob- 
served predominantly when the persistent flux is between 
approximately 2% and 4% of the Eddington limited flux. 
Furthermore, indications were found for the association 
of short recurrence times and the accretion of hydrogen- 
rich material. The shortest recurrence tim e previously 
reported is 5.4 minutes (Linares et al.|[2009l) . 

Different ideas have been put forward to explain this 
rare bursting behavior. As most of the models only re- 
solve the neutron star envelope in the radial direction, it 
is possible that short recurrence time bursts are due to 
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multi-dimensional effects, such as the confinement of ac- 
creted material on different parts of th e surface, possibly 
as the result of a magn eti c field (e.g ., 'Melat os fc Pavng 
l2005HLambet al.ll2009f l. iBoirin et a l. (2007), however, 
found that the different bursts originate from an emit- 
ting area of similar size. Furthermore, the indication of a 
different fuel composition for the bursts with short recur- 
rence times, argues against any scenario where accreted 
material of the same composition burns on different parts 
of the surface. 

The idea of a burning lay er with an unburned layer on 
top has been investigated (jFuiimoto et al.l [19871 ). After 
the first layer flashes, the second layer could be mixed 
down to the depth where a thermonuclear runaway oc- 
curs. Mixing may be driven by rotational hydrodynamic 
instabilities (jFuiimotol Il988f ) or by ins tabilities due to 
a rot ationally induced magnetic field ( Piro fc Bildstenl 
l2007t iKeek ct al. 2009 ) . The mixing processes take place 
on the correct time scale of approximately ten minutes. 
Although this scenario is able to explain many of the 
observed aspects of short recurrence time bursts, it has 
not been reproduced with a multi-zone stellar evolu- 
tion code that inc ludes a full nuclear burning network. 
iTaam et al.l (|1993f ) created models that exhibit 'erratic' 
bursting behavior, reminiscent of short recurrence time 
bursts. Later versions of the employed code, however, 
no longer produce this, most likely because o f the inclu- 
sion o f a more extensive nuclear network ( Wo oslev et al.l 
2004f). As a different explanation for the reignition. 



Boirin et al.l (|2b07) suggested a waiting point in the chain 



of nuclear reactions. There may be a point in the chain 
where a decay reaction with a half life similar to the short 
recurrence times stalls nuclear burning before continuing. 
We use an improved version of the burst catalog com- 
piled by Galloway et al. (2008), that is extended with 
the X-ray burst observ ations of the W FCs on-board Bep- 
poSAX (e.g.. .Cornelisse et al.ll200^ . This is the largest 
collection of X-ray bursts used in any study to date. This 
allows us to study the short recurrence time phenomenon 
in much more detail. 

2. OBSERVATIONS AND ANALYSIS METHODS 

2.1. Nomenclature 

In this paper we name the different kinds of bursts 
using the conventions of IBoirin et al.l ()2007[ ). Bursts 
with recurrence (waiting) times shorter than one hour 
are referred to as short waiting time (SWT) bursts, 
while longer recurrence time bursts are long waiting time 
(LWT) bursts. The one-hour boundary is chosen to dis- 
criminate between the two distinct groups of bursts we 
observe (Sect. l375|) . A burst event is defined as a series 
of bursts, where any two or more subsequent bursts are 
separated by a short waiting time. We refer to an event 
with one, two, three, or four bursts as a single^ double^ 
triple, or quadruple burst, respectively. Furthermore, we 
call a double, triple, or quadruple burst a multiple-hurst 
event. For the bursts within a multiple-burst event, we 
use the terms first hurst and follow-up hurst, where the 
former refers to the LWT burst and the latter to any 
SWT burst in the event. 

2.2. Burst catalog 

Because SWT bursts are rare, we need a large sample 
of bursts to study them well. We use a preliminary ver- 



sion of the Muhi-INstrument Burst ARchive (MINBAR), 
which is a collection of Type I bursts that are observed 
with different X-ray instruments, and that are all ana- 
lyzed in a uniform way. MINBAR is a continuation of the 
effort whi ch started with the RXTE PGA burst catalog 
bv Gallowav et al.l ()2008[ ). combined with the many X- 
ray bursts observed w ith the WFCs on BeppoSAX (e.g., 
iGornelisse et al.ll2003[ ). The catalog will be presented in 
full in a forthcoming paper. Currently it contains in- 
formation on 3402 bursts from 65 sources, among which 
are 1 36 SWT bursts ( MINB AR version 0.4). In compar- 
ison, [cSlQw^ern] (|2008l ) report 1187 bursts from 48 
sources, among which are 84 bursts with a short recur- 
rence time, in that paper defined as < 30 min (MINBAR 
contains 110 SWT bursts using this criterion). 

2.3. Instruments 

The Rossi X-ray Timing Explorer (RXTE) was 
launched on December 30, 1995. One of the in- 
strume nts on-board is the Proportional Counter Array 
(PGA; iJahoda et al.l |2006[ ). consisting of five propor- 
tional counter units (PCUs) which observe in the 1 to 
60keV energy range. The PGA has a large collecting 
area of 8000 cm^ (1600 cm^ for each PCU). The PCUs are 
co-aligned and have a collimator that gives a 1° FWHM 
field of view. As part of the primary mission objective 
of RXTE, the PGA gathered a large amount of exposure 
time on the galactic LMXB population. The observations 
we use have an average duration of 78 min. We use all 
data that was p ublicly available in Ju ly 2008 (compared 
to June 2007 for lGallowav et al.ll2008l) . 

The High-E nergy X-ray Timing Experiment (HEXTE; 
iGruber et al.l [l996) on RXTE consists of two clusters 
of Nal/Csl scintillation detectors that allow for obser- 
vations in the 15 to 250 keV energy range. If available, 
we combine PGA and HEXTE observations to obtain a 
broad band X-ray spectrum, when studying the persis- 
tent fiux of bursting sources. 

The All-Sky Monitor (ASM) on RXTE consists of 
three scanning shadow cameras (SSCs) on a rotating 
beam that image a large part of the sky each satellite 
orbit, in a series of short 90s observations. The SSCs 
observe in the 1.5 to 12keV energy range. 

A few months after RXTE, the BeppoSAX obs ervatory 
was launched in April 1996 (iBoella et al.lll997[ ). and it 
was operational until April 30, 2002. On-board were two 
Wide-Field C amera's (WFC s) which faced in opposite 
directions (Ja ger et al.lll997l ). The WFCs were sensitive 
in the 2 to 28 keV band-pass and each camera imaged at 
any time 40° x 40° of the sky using a coded mask aper- 
ture. Bi-yearly observations of the Galactic Center were 
performed, resulting in a lar ge exposure time for m any of 
the LMXBs in the Galaxy (|in 't Zand et al.ll2004f l. the 
WFG observations have a mean duration of 255 min. We 
use all WFG data. 

For the majority of the bursters, which are located near 
the Galactic Center, we gather a total net exposure time 
of approximately 50 days. 

Both observatories were placed in a low Earth orbit of 
approximately 96 minutes. During the observation of a 
particular source, the source is obscured by the Earth 
for a typical duration of approximately 36 minutes per 
orbit. Furthermore, when the satellites passed through 
the South Atlantic Anomaly (SAA), the detectors were 



X-ray Bursts with Short Recurrence Times 



turned off to prevent damage. This introduces data gaps 
of 13 to 26 minutes long. The precise length of the dif- 
ferent data gaps depends on the latitude of the source 
with respect to the satelhte orbit, which was different 
for BeppoSAX and RXTE. 

2.4. Burst analysis 

We briefly discuss the process of the burst analy- 
sis. The generation of data products for the differ- 
ent instruments is handled identicall y to th e studies by 
iGalloway et al. (2008) and Cornehss eet all (|Ip03). Our 
method of burst analysis is in principle the same as for 
those studies, but extra care was taken to ensure that re- 
sults from different instruments are directly comparable. 
For more d etails of the analysis of PCA and WFC data, 
we ref er to I Galloway et al.l l|2008f l and iCornelisse et al.l 
([2003) , respectively. 

To find burst occurrences, we generate light curves for 
each instrument, for each known bursting source. In the 
light curves we locate all events that rise significantly 
above the background, e.g., exceeding the mean flux level 
by at least four times the standard deviation. These 
events are checked by eye for the characteristic profile of 
a fast rise and an exponential-like decay. 

Next, time- resolved spectroscopy is performed on the 
candidate bursts. We divide the individual bursts in time 
intervals, such that we obtain for each interval a burst 
spectrum with similar and sufficient statistics. The net 
burst spectrum is obtained by subtracting the spectrum 
of the persistent emission observed du ring the entire 
obser vation, excluding the burst (e.g., iKuulkers et al.l 
120021 ). We fit the burst spectra with a black body model, 
taking into account the effects of interstellar absorption 
dMorrison fc McCammojl 119831 ; solar abundances from 
[Anders & Ebihara 1982). From this we flnd a black body 
temperature and radius, as well as the unabsorbed flux. 
By extrapolating the fitted black body model beyond the 
observed energy range, we obtain the bolometric unab- 
sorbed flux. Integrating the flux over the burst yields the 
burst fluence. 

The decay time is determined by flts to the burst light 
curve. We obtain the net burst light curve by subtracting 
the persistent flux, as measured in the entire observation, 
excluding the burst. We start fitting the decay when 
the fiux drops below 90% of the peak flux. This has 
the advantage that we are less sensitive to the Poisson 
noise at the peak or to the effects of radius expansion. 
The decay of many bursts is flt well by two exponentials, 
with two exponential decay times. For some bursts two 
exponentials did not yield a statistically better fit than 
a single exponential. For those bursts we report only a 
single decay time. It should be noted that not for all 
bursts a good fit was obtained, but the 'best' fit still 
provides a qualitative description of the burst decay. 

The determination of the burst decay time is currently 
not done uniformly for the different instruments. For the 
PCA the bolometric flux is used, while for the WFCs the 
flux in counts per second is used. We compared the de- 
cay times obtained for 15 bursts that have been observed 
by both instruments. On average, the decay times for 
the WFC bursts are (22 ± 7)% longer than for the PCA 
bursts. Although the decay time scales are not fully com- 
patible across the different instruments, we can still use 
them to differentiate between short bursts and the longer 



bursts with the rp-process tail. 

2.5. Persistent emission and mass accretion rate 

We obtain the persistent flux for each burst by fitting 
the spectrum from the entire observation, excluding the 
burst. By extending the fitted spectral model (typically a 
comptonized spectru m or a black body + power law; see 
iGallowav et al.|[2008l for details) outside the observed en- 
ergy range, we determine the bolometric correction. The 
uncertainty in the bolometric correction can be as small 
as 10%, when we combine RXTE PCA and HEXTE ob- 
servations to obtain a broadband spectrum. The Bep- 
poSAX WFCs, however, do not have such broad energy 
coverage, and some PCA observations suffer from source 
confusion, which leads to an increased uncertainty in the 
bolometric correction of up to 30%. 

To convert flux to luminosity and fluence to energy, 
we multiply by AttcP, with d the distance to t he so urce. 
We use the distances from IKuulkers et al.l (I2003D for 
globular cluster sources, and from IGallowav et al.l ()2008l ) 
(distances from photosphe ric-radius expans ion bursts ob- 
served with the PCA) and lLiu et"all (2007) for the other 
sources. Most measurements of the distance have an un- 
certainty of the order of 30%. Combined with a 30% error 
in the bolometric correction, this leads to an uncertainty 
in the luminosity and the fluence of up to approximately 
70%. 

The persistent X-ray emission from an LMXB mainly 
originates from the inner part of the accretion disk, and, 
therefore, is a measur e of the mass accretion rate (e.g, 
IGallowav et al.ll2004bf ). We express the mass accretion 
rate M in terms of the Eddington limited accretion 
rate MEdd by equating their ratio to the ratio of the 
persistent luminosity L and the Eddington luminosity 
for hydrogen accreting sources Lsdd = 2 ■ lO'^^ergs"^: 
A//MEdd = L/JjEdd- L/Edd dcpcuds ou the neutron star 
mass and the hydrogen fraction of the accreted mate- 
rial (e.g., lBildstenlll998i ). both of which are not known 
to great precision for most LMXBs . The current obser- 
vatio nal constraints on the mass (jLattimer fc PrakashI 
|2007( ) introduce an uncertainty of several tens of per- 
cents. Furthermore, we assume that the accretion pro- 
cess has an efficiency of 100%. It is possible that part of 
the matter leaves the syste m in a jet, such as observed 
in black- hole binaries (e.g., iFender et al.l 120051 ). In this 
paper, however, we neglect this possibility and assume 
the luminosity is a good measure of the mass accretion 
rate. We also neglect any anisotropy factors that may 
arise from the inclination of the disk with respect to the 
line of sight; because the inclination is ill-constrained for 
most LMXBs, we assume isotropic emission. 

The combined uncertainties are quite large, but this 
only plays a role when we compare different sources. It 
is, however, of no consequence when we compare the 
bursts of any single source. We will still compare dif- 
ferent sources, but one must be careful to keep these 
uncertainties in mind. 

3. RESULTS 

3.1. Source selection 

The MINBAR catalog contains X-ray bursts from 65 
sources, 15 of which exhibit SWT bursts (Table [T]), 
i.e., bursts with a recurrence time shorter than one 
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TABLE 1 
Overview of 15 bursters with short recurrence times. 



Name 



^spm 

(Hz) 



^exposure MINBAR Single Double Triple Qua- Remarks 
(days) bursts druple 



EXO 0748-676 

GS 0836-429 

4U 1323-62 

4U 1608-522 

4U 1636-536 

MXB 1658-298 

4U 1705-44 

XTE J1710-281 

4U 1735-44 

2S 1742-294 

EXO 1745-248 (Tz5)'' 

4U 1746-37 (NGC 6641)^ 

SAX J1747.0-2853 

Aql X-1 

Cyg X-2f 



552=^ 


83.1 


269 




30.9 


17 




53.9 


41 


620'' 


52.4 


67 


581== 


60.6 


241 


567^* 


49.8 


27 




60.1 


125 




47.9 


18 




54.7 


50 




56.3 


269 




48.3 


24 




55.1 


31 




71.5 


63 


549° 


34.3 


60 




55.7 


55 



251 


9 







Triples: Boirin et al. f2007): ApDendix 


17 









Double: Aoki et al. fl992) 


33 


4 








60 


2 


1 






212 


11 


1 


1 




25 


1 








108 


5 


1 


1 




16 


1 








38 


6 








240 


10 


3 






15 


1 


1 


1 


Type II bursts? CGalloway et al. 2008) 


15 


8 






Two sources? ("Galloway et al. 2004a') 


57 


3 








49 


4 


1 






45 


5 









=" iGallowav et al.l Il2009f) 
''[Muno et aU f2Q02 ) 
" [Strohmayer e t aTTllT99l) 
=* [Wiinands et atl S20mt) 
" [Zhang et al.l 1 11998) 
Source excluded from our analysis. 



See Sect. l3.1l for details. 



hour. We exclude the Rapid Burster, which is known 
to exhibit Type I I bursts, wh i ch ar e not of thermonu- 
clear origin (e.g.. iLewin et ai1ll993( ). Interestingly, no 
multiple-burst events are detected from candidate and 
confirmed ultra-compact binaries (UCXBs). The com- 
panion star in a UCXB is thought to be an evolved 
star, donating hydroge n-poor matter to the neutron star 
(|in 't Zand etaij 120050 . For confirmed UCXBs the bi- 
nary period has been measured, while candidates are 
identified by tentative measurements of the period, by 
a low optical to X-ray flux, or by stable mass transfer at 
rates below 1% of the Eddington limit ed rate. We em- 
ploy the list of (candidate) UCXBs from lin 't Zand et al.l 
|2(307), which omits candidates proposed on the ba- 
sis of their X-ray spectrum, which is likely a less reli- 
able method. From these sources we find 229 bursts, 
none of which are SWT bursts. The frequent burster 
4U 1728-34 (GX 354-0) is suspe cted of being a UCXB 
based on its bursting behavior ([Galloway et al.l [2008D . 
543 bursts of this source are present in MINBAR, but 
no multiple-burst events. As this strongly suggests that 
SWT bursts are limited to hydrogen-rich accretors, we 
exclude in our studies t he list of (candidate) UCXBs from 
lin 't Zand et al.l (|2007D . with the addition of 4U 1728-34. 

GX 17+2 and Cyg X-2 exhibit b ursts at accretion rate s 
close to the Eddington limit (e.g.. lKuulkers et al.ll200^ . 
while most LMXBs do not sho w bur sts above approxi- 
mately 10% of Edding ton (e.g.. lvain^aradiis et al. 1988.; 
ICornelisse et al.l [20031 ). Of these two sources, Cyg X-2 
exhibits SWT bursts. Due to the high accretion rate, 
however, the amount of matter accreted in between the 
bursts is enough to account f or the amount of fuel bu rned 
in the bursts. Furthermore, I Galloway et al.l (|2008D find 
indications that (some of) the bursts of GX 17+2 and 
Cyg X-2 could be Type II bursts. For these reasons we 
exclude these two sources as well. 

The bursts from EXO 1745-248, located in the globular 
cluster Terzan 5, exhibit only weak evidence for cooling, 
which means that their thermonuclear origin is not firmly 



established, and that they possibly are Type II bursts. 
Another sign that the bursting behavior is anomalous, is 
the fact that most sources in Table [1] exhibit more double 
bursts than triples and quadruples, while EXO 1745-248 
has one of each. We exclude the source in our studies. 

4U 1746-37, located in the globular clu ster NGC 6641 
exhibit s both faint and bright bursts. I Galloway et al.l 
(l2004aD found in PCA observations that the faint bursts 
occur at very regular intervals, unaffected by the occur- 
rence of bright bursts, and vice versa. This lead to the 
speculation that the faint bursts originate from a dif- 
ferent LMXB that is also located in NGC 6641. We, 
therefore, exclude the bursts from this source. 

After excluding the mentioned sources, we consider 44 
hydrogen-rich accretors from which we observe 2274 sin- 
gle, 56 double, 7 triple, and 2 quadruple events. 

3.2. Source confusion 

While the WFCs are imaging instruments, the PCA is 
not. The PCA's collimators restrict the field of view 
to 1° FWHM. Especially in crowded regions, such as 
near the Galactic Center, multiple X-ray sources may 
be in the field of view. For the sources in Table [T] 
we check the RXTE ASM light curves for any nearby 
bright X-ray sources that were active at the time of PCA 
burst observations. This is the case for 2S 1742-294 and 
SAX J1747. 0-2853. Those bursts, for which we cannot 
reliably measure the persistent flux, are excluded from 
our studies. 

The problem is smah for SAX J1747.0-2853, because 
most of its bursts were observed with the WFCs, and we 
have a reliable measurement of the persistent flux for all 
three SWT bursts. For 2S 1742-294 the problem of source 
confusion plays a role for 61 bursts, including 12 out of 
16 SWT bursts. Most of these bursts, including all SWT 
bursts, occur in the time interval MJD 52175-52195. The 
persistent flux as measured with the WFCs for bursts 
from that period varies by less than 10%. We assign the 
mean WFC persistent flux to those PCA bursts. 
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Fig. 1. — Histogram of simulated recurrence times irecur • The in- 
trinsic distribution is from a Monte Carlo simulation of 10® bursts 
with SWT and LWT recurrence times as well as SWT fraction from 
EXO 0748-676 (Boirin et al. 2007). The observed distribution in- 
cludes the effect of data gaps due to Earth occultation and the 
South- Atlantic Anomaly. AH bursts with trocur < 60 min. (dotted 
line) are considered SWT bursts. 



3.3. Data gaps and recurrence time 

We define the recurrence time as the time since the pre- 
vious burst from the same source in the catalog. Due to 
the frequent data gaps (Sect. [273)) . we must keep in mind 
that these are upper Umits. Performing Monte Carlo 
simulations, we investigate the effect of the data gaps on 
the number of observed SWT and LWT bursts. 

We generate a series of 10^ burst occurrence times and 
check which bursts fall in data gaps. We use EXO 0748- 
676 as a template: we generate LWT and SWT bursts 
with an SWT fraction of 30%; we use the mean LWT 
and SWT re currence times 3.0 hours and 12.7 minutes, 
respectively (|Boirin et al.l 120071 ): we position the bursts 
in time following a Gaussian distribution around the 
LWT or SWT t,.ecur, with a width of 16% of either t^cur 
(mean variabili ty in pe rsistent flux of a series of per- 
sistent sources; iKeek et al. 2006) to model variations in 
the mass accretion rate. To check which of these bursts 
would be observed in the presence of data gaps, we as- 
sume a 96 min. satellite orbit, containing a 36 min. data 
gap due to Earth occultation. The presence and duration 
of data gaps due to the South- Atlantic Anomaly (SAA) 
depends on the position of the source and the satellite at 
the time of the observation. We model SAA data gaps by 
placing a 20 min. gap at a random phase of each orbit. 

Only 49% of the generated bursts is 'observed', the rest 
coincide with a data gap. When a burst is missed, the 
recurrence time of the next burst, as seen from the pre- 
viously detected burst, is incorrectly found to be longer: 
in the distribution of observed recurrence times a long 
tail of bursts is present at longer trccur than present in 
the original distribution (Fig. [1]). There are also bursts 
in between the Gaussian peaks of SWT and LWT bursts, 
but their number is less than 10% of the total num- 
ber of SWT bursts. Because most of these bursts have 
irecur < 60 miu., wc Still considcr them SWT bursts. 

The fraction of bursts that are SWT bursts, is reduced 
when one or more bursts in a multiple-burst event occur 
during a data gap. While we start our simulations with 
an SWT fraction of 30%, the observed distribution has a 
fraction of 20%. We repeat the simulations with different 
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Fig. 2. — SWT fraction as a function of the SWT recurrence 
time trccur ■ The solid line allows for variation of trocur following a 
Gaussian distribution, while the dashed line does not. The latter 
drops to at irecur = 60 min. , which is the largest SWT trccur we 
consider. 
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Fig. 3. — SWT fraction as a function of the duration of the 
Earth-occulation data gaps. The SWT fraction drops to when 
the data gaps approach the 96 min. duration of the satellite orbit. 

values of the SWT trccur- The total fraction of detected 
bursts remains the same, but the SWT fraction drops 
from 30% to 10% for increasing trccun until trccur exceeds 
the duration of the Earth occultation data gap (Fig. [5]). 

Repeating the simulations with different durations of 
the Earth-occultation data gap, the SWT fraction drops 
by only a few percent for longer data gaps, until the 
fraction quickly goes to 0, when an SWT recurrence time 
no longer fits in the observed part of the orbit (Fig. [3]). 

A similar decrease of the SWT fraction is expected if 
the duration of the observation is less than the SWT re- 
currence time. While the average exposure times of both 
the PCA and the WFCs exceed one hour, a substantial 
part of the observations had a shorter recurrence time: 
19% of the WFC exposures and 66% of the PCA expo- 
sures were shorter than one hour. 

3.4. Detection limits 

The instruments we use have different detection limits. 
The WFCs had a substantially larger field of view than 
the PCA, which results in a higher background level. 
Furthermore, the PCA has a 46 times larger collecting 
area than each WFC. Consequently, we are able to find 
fainter bursts in PCA data than in WFC data (Fig. g]). 
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Fig. 4. — Histogram of observed peak burst flux -Fp^ak for the 
WFCs and the PCA. Additionally we show the distributions for 
the short recurrence time bursts. 

This is especially important for the SWT bursts, as they 
have been fo und to be o n average fainter than the LWT 
bursts (e.g., iBoirin eFlil. 20071. In the PCA data we 
find SWT bursts with peak flux Fpoak as low as 1.6 • 
10~^°ergcin~^s~^, while in the WFC data the faintest 
SWT burst has Fpoak = 3.6 ■ lO-^ergcm-^ s'^ As a 
result, we find many more SWT bursts with the PCA: 
for the PCA we find 76 SWT bursts out of a total of 910 
bursts, and for the WFCs we find 14 SWT bursts out of 
1560 bursts. 

3.5. Recurrence times 

We plot for all bursts the recurrence time irecur as 
a function of the persistent luminosity Lpers (Fig. [S]). 
While most bursts have a recurrence time of at least sev- 
eral hours, there is also a group of SWT bursts with 
^recur < 1 hour. There is an intrinsic spread in trccur due 
to, for example, variations in the mass accretion rate, or 
variations in the temperature in the neutron star enve- 
lope. We investigate whether the short recurrence times 
can be explained as the tail of the distribution of the 
long recurrence times. For this distribution we assume 
a Gaussian, even though it is not certain whether this 
is correct far from the mean. The data gaps modify 
the observed distribution, especially towards longer irccur 
(Sect. l575|) . The leading part of the Gaussian, however, is 
not modified substantially, apart from the overall lower 
number of observed bursts due to the lower net exposure 
time (Fig. [T]). We fit a Gaussian to the distribution of 
recurrence times between 1 and 3 hours with the cen- 
ter fixed at 3 hours. Extrapolating the best fit towards 
shorter irccur, we predict 5.6 SWT bursts. The Poisson 
probability for the observed number of SWT bursts of 
76, is negligibly small {P < 10""'^°). Therefore, the short 
recurrence times follow a separate distribution. 

There is a separation between bursts with short (< 
0.5 hour) and long (> 1 hour) recurrence times. Above 
-^pors ^ 6 • lO^^ergs^^, however, there are some bursts 
that have recurrence times of 30 to 60 minutes. Com- 
paring the recurrence time distributions of SWT bursts 
(irecur < 1 hour) at persistent luminosities below and 
above Lpers — 6 • lO^^ergs"-^, a KS-test yields P = 0.16, 
which means we can exclude at 84% that both distri- 
butions are the same. This is not a strong constraint 
constraint, mainly due to the small number of bursts 



with Lpcrs < 6 • 10^^ ergs '^. It is, however, consistent 
with the behavior of EXO 0748-676 during EXOSAT, 
XMM-Newton, and Chandra observations, that resulted 
in relatively large data sets, where SWT bursts with re- 
currence times exceeding 30 minutes only occurred when 
the persistent flux was lar ger (jGold 1968 : Gottwald et al.l 
ll987bl:[Boinn et al.ll2007l : Appendix). 

There are SWT bursts at all values of Lpers where 
LWT bursts are observed, with the possible exception of 
Lpeis ^ 3- lO^^ergs"^, although this may be a statistical 
effect due to the lower number of bursts. 

The a-parameter is defined as the ratio of the persis- 
tent fiuence between bursts and the burst fiuence. As- 
suming a burst fiuence of 2.2 ■ 10'^^ erg — the average 
fiuence of MINBAR single bursts from hydrogen-rich ac- 
cretors — we draw lines of constant a; a ~ 40 is a typical 
value for many burs ters, while a o:^ 1 000 is observed for 
superbursters (e.g.. Hn 't Zand et al.li2003 ). The effec- 
tive a- value for the SWT bursts is far below a = 40, 
which is the expected value for thermonuclear ignition 
of mixed hydrogen/helium fuel, highlighting the require- 
ment for ignition of unburned fuel left-over from the pre- 
vious burst. 

We find the shortest recurrence time reported so-faiQ: 
3.8 minutes for a double burst from 4U 1705-44, detected 
with the WFCs at MJD 51233.89. 

3.6. Accretion rate dependence 

We use the persistent luminosity Lpcrs as a measure 
of the mass accretion rate (Sect. 12. 5|) . A Kolmogorov- 
Smirnov (KS) test finds that the distributions of all SWT 
and LWT as a function oiLpas are compatible (P = 0.10; 
Fig. [5]): there are multiple-burst events at all mass ac- 
cretion rates where normal bursts occur. This does not 
hold, however, for all individual sources. We investigate 
a few frequent bursters more closely. For EXO 0748-676 
the distributions of single and multiple bursts are com- 
patible (P = 0.73); 4U 1636-53 and 2S 1742-294 exhibit 
multiple bursts only in a small Lpers interval, while sin- 
gle bursts occur in a wider range of Lpcrs- The small 
Lpcrs intervals for these two sources do not seem to coin- 
cide. The uncertainty in the luminosity, however, is large 
(Sect. 12. 5p . so the intervals may still be consistent. 

We find SWT bursts for only 15 out of 44 hydrogen-rich 
accretors. For most of the sources we can attribute the 
lack of SWT bursts to the low number of bursts detected 
per source. There are, however, two bursters, KS 1731- 
260 and GS 1826-24, from which we have detected over 
300 LWT bursts per source, but no SWT bursts. 

The position in so-called color-color diagrams, Sz, 
is regarded as a tracer of t he mass accretion rate 
(jHasinger fc van der Klij|1989r ). We compared the dis- 
tribution of Sz for LWT and SWT bursts observed with 
the PCA from eight sources for which Sz is well defined: 
4U 1608-522, 4U 1636-536, 4U 1702-429, 4U 1705-44, 
4U 1728-34, KS 1731-260, Aql X-1, and XTE J2123- 
058, omitting 4U 17 46-37 as explained in Sect. 13.11 
(|Gallowavet al.ll2008[ ). While LWT bursts have associ- 
ated Sz values of up to 2.8, SWT bursts all have 5z < 2. 



* IWiinands et al.l 1)20021 ) found a pair of candidate bursts from 
MXB 1659-298 only ~ 50 s apart. The first flare, however, is rela- 
tively weak and lacks the cooling characteristic for Type I bursts. 
The second flare has been identified as 7-ray burst GRB 990419C 
(HEASARC IPNGRB catalog). 
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Fig. 5. — Observed recurrence time trccur as a function of the persistent luminosity Lpors for 2415 bursts from 44 LMXBs. Due to the 
presence of data gaps, values of trccur exceeding one hour are upper limits to the real recurrence time. We show bursts with trccur < 10'' hour. 
The lines represent constant values of the average a-parameter as indicated (see text). Note that the crosses indicate only the positions of 
the data points, not the uncertainties. The shortest recurrence time is 3.8 minutes. 



Therefore, we find that SWT bursts are restricted to the 
so-caUed island state, while LWT bursts also occur in the 
'banana' branch. A KS-test yields P ~ 10"^, confirm- 
ing that the Si distributions for LWT and SWT bursts 
are different. Most of the SWT bursts from the fre- 
quent burster 4U 1636-53 occurred with 1.5 "^ Sz < 2.0. 
KS 1731-260 exhibits only a few LWT bursts in that 
range. This suggests that SWT bursts occur mainly in a 
small Sz interval, and that the lack of SWT bursts from 
the latter source is caused by the low number of observed 
bursts in that range. Note that we observe SWT bursts 
with Sz as low as 0.8, so the interval is not the same for 
all sources. 

3.7. Frequency of short vs. long recurrence times 

We consider 2415 bursts from hydrogen-accreting 
sources. 76 have a short recurrence time: the overall 
SWT fraction is (3.1 ± 0.4)%. The l-cr uncertainty is 
derived from the Poisson uncertainties in the number of 
(SWT) bursts. As a function of persistent luminosity 
there is some variation in the SWT fraction, but there 
is no clear trend (Fig. [7]). The weighted mean of the 
SWT fraction in the range of Lpers where SWT bursts 
are observed is (2.3 ± 0.3)% for all hydrogen-rich accre- 
tors, (13±4)% for 4U 1636-536, (16±4)% for 2S 1742-294 
and (6 ± 2)% for EXO 0748-676. 

We compare these fractions to XMM-Newton and 
Chandra observations of EXO 0748-676. In 2003 the 
XMM-Newton EPIC PN observed double and triple 
bursts from EXO 0748-676 with an SWT fraction of 



(32 ± 7)% (|Boirin et al.|[200l . IHoman et all ((200l find 
in XMM-Newton EPIC PN and MOS observations from 
2000 and 2001 4 single and 4 double bursts, which results 
in an SWT fraction of (33 ± 19)%. Chandra ACIS-S ob- 
servations from 2001 and 2003 exhibit 41 bursts with an 
SWT fraction of (27±9)% (see Appendix). The weighted 
mean of these rates is (30 ± 5)%. An important differ- 
ence in the data of on the one hand XMM-Newton and 
Chandra, and on the other hand RXTE and BeppoSAX, 
is the frequent data gaps in observations of the latter 
observatories. From Monte Carlo simulations we found 
that data gaps reduce an SWT fraction of 30% to 20% 
(Sect. 133]). This is significantly higher than the (6 ±2)% 
we obtain from MINBAR. We repeat the simulations for 
4U 1636-536 and 2S 1742-294, using the mean SWT re- 
currence times from MINBAR: 16.2 min. and 20.5 min., 
respectively. The obtained SWT fractions, respectively 
18% and 16%, are consistent within 1.3a with the frac- 
tions from MINBAR. 

The discrepancy for EXO 0748-676 may arise because 
of the fact that the WFCs are less sensitive to fainter 
bursts than the PCA or the instruments on XMM- 
Newton and Chandra. Taking only the PCA bursts 
into account, we find SWT fractions of (11 ± 6)% for 
EXO 0748-676, (13 ±4)% for 4U 1636-53, and (23 ±6)% 
for 2S 1742-294, aU of which are within 1.5a from the 
fractions found from the Monte Carlo simulations. 

Two frequent bursters exhibit no SWT bursts: 
KS 1731-26 and GS 1826-24. Combining the number 
of LWT bursts observed from these sources, we derive an 
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Fig. 6. — Histogram of bolometric persistent luminosity Lpors of LWT and SWT bursts, for all hydrogen-rich accretors (top-left) and for 
five individual sources. There is a dip in the distribution for all sources around Lpcrs — 1.7 ■ 10^7 ergs~^. This is due to similar dips in the 
distributions of 4U 1636-53 and 2S 1742-294, as well as peaks in the distributions of KS 1731-26 and GS 1826-24 at slightly lower Lpcrs- 



upper limit to the SWT fraction of 0.14%. This is over 
20 times smaller than the SWT fraction we find for all 
hydrogen-accreting sources combined. 

3.8. Temperature and energetics 

From time-resolved spectral analysis of the bursts, we 
obtain the black-body temperature and the burst en- 
ergetics. The peak temperature of the first bursts in 
multiple-burst events is on average higher than the peak 
temperature of the SWT bursts (Fig. [8]). A KS test shows 
that the temperature distributions for single bursts and 
first bursts are not compatible (P ~ 10~^^). 

We compare the peak luminosity and the fluence of 
single bursts to those in multiple-burst events (Fig. [9l 



Fig. fTO]) . The SWT bursts have on average a lower peak 
luminosity and a lower fluence. The energetics of the first 
bursts does not follow the same distributions as the single 
bursts {P < 10~^). By adding the fluence of all bursts 
in a multiple-burst event, we calculate the total event 
fluence (Fig. [TT|) . Multiple-burst events are on average 
more energetic than single burst events, but do not have 
a higher fluence than the most energetic single bursts. 

Summarizing, SWT bursts are on average weaker and 
cooler than LWT bursts, but the combined fluence in 
multiple-burst events is on average 8% higher than the 
fluence of single-burst events. 

3.9. Decay time scale 
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Fig. 7. — The fraction of bursts that have a short recurrence time as a function of the persistent luminosity Lpers, for all hydrogen-rich 
accretors (top-left) and for three individual sources. For each plot Lpcrs is divided into 30 logarithmically spaced bins, and for each bin the 
SWT fraction is calculated from the bursts in the corresponding plots in Fig. [6] For EXO 0748-676 one SWT burst is observed at higher 



^pors, which results in a data point with large uncertainty. 
103 



10^ 



CD 



lOi 



10" 



T 



1 r 



-r 



T 



-r 



T 




;] Single 
3 First 



r.v.'.i SWT 



f;-^ . -\ 



n 



J I I I L 



lO'^ 



10^ 



a; 

X2 

a 



10^ 



10" 




0.0 0.5 1.0 



1.5 2.0 2.5 3.0 
Peak kr (keV) 



3.5 4.0 4.5 



1036 



1039 



Fig. 8. — Histogram of peak black-body temperature of single 
bursts, the first bursts in multiple events, and the SWT bursts. 

The two component exponential provides a good fit 
to many bursts. If the two-component decay does not 
provide a significantly better fit than the one-component, 
we use the latter. We compare the longest decay time 
scales of all bursts in Fig. [T^l On average, the SWT 
bursts decay faster than the other bursts. KS tests show 
that all distributions are incompatible (P < 10^^). 

3.10. Neutron stars with SWT bursts spin fast 

The spin frequency i^spin of an accreting neutron star 
is measured in observations of accretion-powered pulsa- 
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Fig. 9. — Histogram of bolometric peak luminosity of single 
bursts, the first bursts in multiple events, and the SWT bursts. 



tions or burst oscillations. Both mechanisms are thought 
to arise from hotter and, hence, brighter spots on the sur- 
face rotating in and out of view. Currently :^spin is known 
for 2 5 accreting neut ron stars, 19 of which are bursters 
(e.g., I Gallowavl 120081 ) . Among them are 5 sources that 
exhibit SWT bursts (Table [Ij . They are concentrated 
towards the high-frequency part of the distribution for 
all bursting sources (Fig. [T5)) . Since we only have five 
multiple-bursting sources with a known spin frequency, 
we cannot exclude that this bias is the result of the small 
sample. 
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Fig. 10. — Histogram of bolometric fluence of single bursts, the 
first bursts in multiple events, and the SWT bursts. 
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Fig. 12. — Histogram of the exponential decay time for single 
bursts of hydrogen-rich accretors, as well as the flrst and SWT 
bursts in multiple events. 



There could be a selection effect: sources with a higher 
mass accretion rate accrete more angular momentum, 
causing them to spin up faster. Furthermore, their burst 
rate is higher, making it easier to detect a rare multiple- 
burst event. We check the MINBAR catalog for the num- 
ber of bursts from these sources as a function of the spin 
frequency (Fig. I13p . There is roughly an equal number 
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Fig. 13. — Histogram of the neutron star spin Vgpin determined 
from X-ray observations, for all LMXBs with known spin, for the 
known bursters and for the sources that exhibit multiple-bursts 
(Table [TJ. The latter are concentrated at the high end of the 
distribution of known spins. The dotted line indicates the number 
of bursts in the MINBAR catalog for the bursters in each bin (right- 
hand axis). 
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Fig. 14. — Quadruple burst from 4U 1636-53 as observed with 
the RXTE PCA on M JD 52286. Light curve at 2 s time resolution 
{top) and zoomed in on each burst at 1 s time resolution {bottom). 
For each burst we indicate the (longest) exponential decay time r. 

of bursts observed at higher and at lower i^spin, so the 
selection effect is not present. 

There are a few hydrogen-accreting sources of which 
we detected a large number of bursts, but no SWT 
bursts (Sect. 13. 6p . One of these sources, KS 1731-260, 
is known to have a large spin frequency of I'spin = 524 Hz 
IjSmith et al.lll997| ). 

3.11. Quadruple hurst from ^U 1636-53 

To illustrate the properties of LWT and SWT bursts, of 
which we have shown the distributions using a large num- 
ber of bursts from the catalog, we consider the quadru- 
ple burst from 4U 1636-53 (Fig. [T4|) . The four bursts 
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occurred within 54 minutes, and the time between the 
burst onsets is 18.2 min, 17.9 min, and 16.8 min, respec- 
tively. As indicated in Fig. [Ml the first burst has by 
far the longest decay time r. The first burst has a 
peak flux and fluence that is over seven times larger 
than any of the three SWT bursts. The peak black- 
body temperature (kT) of the four bursts is, respectively, 
(1.67±0.02)keV, (1.06±0.06)keV, (1.51±0.06)keV, and 
(1.50 ± 0.04) keV: again the first burst has the highest 
value. The combined net burst fluence of the quadruple 
event is (2.93 ± 0.05) lO^^erg, which is 35% higher than 
the average fluence of single bursts from this source, and 
27% higher than the fluence of the first burst from the 
quadruplet. This means that at least 27% of the available 
fuel did not burn in the initial burst. 

4. DISCUSSION 

We study the short recurrence time behavior in a large 
sample of bursts from multiple sources. We use a prelim- 
inary version of the MINBAR burst catalog, containing 
3387 Type I X-ray bursts from 65 sources. 15 sources 
exhibit bursts with recurrence times less then one hour: 
SWT bursts. The short recurrence times do not allow for 
the accretion of the hydrogen and helium that is burned 
during the burst, which means that it must have been 
accreted before the previous burst. For example, the 
burst fluences from the quadruple event we observe from 
4U 1636-53 indicate that at least 27% of the accreted 
fuel did not burn in the first burst. This is in contra- 
diction with current one-dimensional multi-zone models, 
which predict that durin g a flash over 90% o f the avail- 
able fuel is burned (e.g.. IWooslev et al.ll2004l ). The aim 
of this study is to provide a comprehensive observational 
assessment of the SWT behavior. 

4.1. Temperature and energetics 

We perform time resolved spectroscopy of the bursts, 
and flnd that SWT bursts are on average less bright and 
reach a lower black-body temperature at the peak than 
the LWT bursts. The fluence of SWT bursts is on aver- 
age lower, but the combined bursts in a multiple-burst 
event are as energetic as the most energetic single bursts. 
This is in agreement with the XMM- Newton observa tions 
of EXO 0748-676 analyzed by Boi rin et all ()2007[) . In 
contrast to that investigation, however, we find that the 
distributions of these quantities for single and for first 
bursts of multiple-burst events are not compatible ac- 
cording to Kolniogorov-Sniirnov tests. It is possible that 
our improved statistics allows us to see a disparity that 
previously went unnoticed. 

4.2. Multiple surface regions vs. multiple layers 

Several ideas have been put forward to explain SWT 
bursts in an attempt to answer the two main questions: 
how to preserve fuel during a burst for the next burst, 
and how to reignite this fuel on a time scale of approxi- 
mately ten minutes. Subsequent bursts with short recur- 
rence time may take place in different regions of the neu- 
tron star surface. The accreted matter may be confined 
by a magnetic dipole field to the poles of the neutron 
star, or perhaps the burning front of a burst is stalled 
at the equator if the inflow of matter from the accretion 
disk is particularly strong there. This would provide an 



explanation for double bursts, but the triple and quadru- 
ple bursts that we observe would require a more compli- 
cated configu ration of the magnetic field. Furthermore, 
IBoirin et al] (pOOl find for EXO 0748-676 that there is 
no evidence for a difference in the X-ray emitting region 
during the different bursts of multiple-burst events. Also, 
they find indications that SWT bursts burn a fuel mix- 
ture with a lower hydrogen content, which would not be 
the case if different regions of pristine accreted material 
are burned. We, therefore, favor the scenario where SWT 
bursts take place in d ifferent layers on top of each other 
(jFuiimoto et al.iri987f) . For this to work, the thermonu- 
clear burning during the flash must be halted before the 
hydrogen and helium is depleted. 

4.3. No SWT bursts from UCXBs 

We observe no SWT bursts from any of the 16 (can- 
didate) ultra-compact sources from which in total 229 
bursts have been observed with the PCA and WFCs. For 
this reason we excluded these sources from our analyses. 
For the hydrogen-accreting sources we found an fraction 
of 3.1% of all bursts are SWT bursts. If we assume 
the same SWT fraction for the most frequently burst- 
ing conflrmed UCXB, 4U 1820-303, we expect 1.7 SWT 
bursts out of the 54 bursts observed by the PCA and 
WFCs. The Poisson probability of detecting no SWT 
bursts when expecting 1.7, is 0.18. Taking into account 
the bursts from all confirmed and candidate UCXBs, 
we expect 7 SWT bursts, and the probability of a non- 
detection is less then 10"'^. 

Based on its burstin g behavior, 4 U 1728-34 is a sus- 
pected UCXB (jGallowav et al.ll2008D . If we include the 
543 observed by the PCA and WFCs from this source, 
the expected number of SWT bursts is 23, and the prob- 
ability of detecting none is 10^^°. 

4.4. Thermonuclear burning processes 

No SWT bursts are observed from UCXBs, but only 
from sources that are thought to accrete hydrogen-rich 
matter. This suggests that the nuclear burning pro- 
cesses involving hydrogen, i.e., the hot CNO cycle, the 
ap-process, and the rp-process, are important for cre- 
ating SWT bursts. The rp-process is a series of proton 
captures and /3-decays, that creates heavy isotopes with 
mass numbers up to approximately 100 (jSchatz et al.l 
120011 ). In this reaction chain there might be a nuclear 
waiting point with the correct time scale to interrupt and 
reignite the thermonuclear burning, for example the time 
scale for spontaneous /3-decay of an isotope. We find, 
however, a broad distribution of short recurrence times 
ticcuT (Fig- El J which argues against a single waiting point 
in the reaction chain (see also IBoirin et al.ll2007[) . Note 
that we do not detect SWT bursts from all hydrogen- 
rich accretors, the two frequent bursters KS 1731-26 
and GS 1826-24 being the best examples. This means 
that merely accreting hydrogen is not enough to produce 
SWT bursts. 

The decay profile of an X-ray burst is shaped by two 
processes. First there is radiative cooling on a ther- 
mal time scale. For normal bursts, which ignite at a 
typical column depth of y ~ lO^gcm'^, this time scale 
is Tthcrm — 10 s. A sccoud process, that slows down 
the decay, is prolonged thermonuclear burning through 
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the rp-process, wh ich lasts up to approxi mately 100 s 



Schatz et alj mmh . For EXO 0748-676 iBoirin et all 



20071 ) found SWT bursts to lack the second slower de- 



cay component, while the first bursts clearly exhibit a 
two-component exponential decay. We confirm that this 
holds true for the other sources with SWT bursts as well. 
This supports the conclusion by IBoirin et al.l (|2007ll that 
follow-up bursts must occur in a layer with a significantly 
r educed hydrogen c onten t . 

lin 't Zand et al.l ([2003) found observations where the 
burst decay can be followed for several thousands of sec- 
onds. They explain the long tail as due to the cooling of 
a deeper layer below the bursting layer, that is heated by 
the burst. One of these long tails is detected for the first 
burst in a triple event of EXO 0748-676. The tail con- 
tinues to decay uninterrupted while the second and third 
bursts occur. This is consistent with the idea that SWT 
bursts occur in a layer above the ignition depth where 
the first burst occurs. Taking into account that SWT 
bursts are less energetic, the deeper layer would not be 
heated substantially by the SWT bursts and continues 
to cool. 

4.5. Mass accretion rate dependence 
SWT bursts are observed over the entire range of mass 
accretion rates M where LWT bursts are observed. For 
some individual sources, however, the SWT bursts occur 
in a smaller M interval than LWT bursts. This is sup- 
ported by the position in the color-color diagram where 
SWT bursts are observed (Sect. 13. 6| Gal loway et al.l 
120081 ). Other frequent bursters exhibit no SWT bursts at 
all, even though the ranges of M we observe from sources 
with and without SWT bursts overlap. Because of the 
large uncertainty in converting flux to accretion rate, the 
precise overlap is uncertain. 

At low accretion rates SWT recurrence times irccur 
are mostly restricted to 3.8 min < ^rccur ^ 40 min. 
Above approximately 0.05 MEdd, where the Edding- 
ton limited mass accretion rate Msdd corresponds to 
a persistent luminosity of Lsdd = 2 • lO^^ergs"^ for 
hydrogen-accreting sources, recurrence times occur also 

in the range 40 min < iiocur ^ 60 min. At 0.05 MEdd a 
transition between tw o bursts regimes is pr edicted by 
iFuiimoto etHI (|1981D (see also lBildstenlll998l ). For lower 
accretion rates all accreted hydrogen burns in a stable 
manner, and the burst ignites in a hydrogen-poor layer. 
At higher rates there is no time to burn all hydrogen, 
and the burst ignites in a layer containing a substantial 
fraction of both hydrogen and helium. It may be that the 
latter regime allows for short recurrence times as long as 
an hour, while the former regime does not. 

4.6. Rotation and mixing 

The spin frequency I'spin is not known for most accret- 
ing neutron stars, as it requires the observation of X-ray 
pulsations or burst oscillations (e.g. . IGalloway|200§ l . For 
five sources with short recurrence times i/gpin is known: 
all five are fast spinning neutron stars with t'spin ^ 500 Hz 
(Table [IJ . The fast rotation could be required for the 
occurrence of multiple-burst events. It induces rota- 
tional instabilities, for example shear instabilities (e.g. 
iFuiimotO : ,1988 ). and instabilities due to a rotationally 



induced magnetic field (|Spruitl[2002D . that mix the neu- 
tron star envelope on a time scale of approxima tely ten 
minutes (jPiro fc Bildstenl 120071: Ikeek et al.ll2009l) . If the 
thermonuclear burning during a flash is halted before it 
reaches higher layers, the hydrogen and helium in those 
layers will be mixed down. On a ten minute time scale 
it reaches the depth where temperature and density are 
sufficiently high to create the thermonuclear runaway for 
the next burst. 

At accretion rates higher than 0.05 AfEdd, we find 
bursts with short recurrence times as long as an hour. 
We mentioned that this may be related to the tran- 
sition to a different burst regime. The time scale for 
rotational mixing depends strongly on the thermal and 
composition al profile of the neutron star envelope (e.g., 
IHeger et"~ai . 2000 for the case of massive stars), which 
may vary for different burning regimes or even differ- 
ent bursts. This could provide an explanation for the 
spread in the observed recurrence times. Further theo- 
retical study is necessary to better understand this. 

The occurrence of multiple-burst events in sources with 
high rotation rates has consequences for the hypothesis 
that strong magnetic fields at the neutron star surface 
contain accreted matter at the poles. This would ex- 
plain multiple bursts as caused by the burning of dif- 
ferent magnetically-confined patches at the poles. The 
presence of a strong magnetic field, however, would al- 
low for the transportation of angular momentum away 
from the neutron star, causing it to spin slower. The ob- 
servations of short recurrence time bursts preferentially 
at high fspin seems in contradiction with this, which dis- 
favors the magnetic-confinement scenario. 

KS 1731-260 spins at a high frequency of 524 Hz, 
accreted hydrogen-rich material, and exhibited many 
bursts (369 in MINB AR) . No short recurrence times were 
observed. Therefore, while a high rotation rate may sup- 
port the occurrence of SWT bursts, it is not possible 
to discriminate between sources with and without SWT 
bursts based on this property alone. Possibly a combina- 
tion of fast rotation and a mass accretion rate within a 
certain range (see previous section) are required for short 
recurrence times. 

4.7. Frequency of SWT bursts 

We investigate the SWT fraction: the number of bursts 
that have a short recurrence time with respect to the to- 
tal number of observed bursts. We determine this frac- 
tion, at the persistent luminosities where SWT bursts are 
observed, for three frequent bursters with SWT bursts: 
EXO 0748-676, 4U 1636-53, and 2S 1742-294. XMM- 
Newton and Chandra observations of EXO 0748-676 find 
an SWT fraction of (30 ± 5)%. Our burst sample is ob- 
tained from RXTE PGA and BeppoSAX WFG observa- 
tions, which contain data gaps due to Earth occultations 
and due to the South-Atlantic Anomaly. Monte-Carlo 
simulations show that the data gaps reduce an SWT 
fraction of 30% to 20%. An additional problem is the 
fact that the WFCs are less sensitive to fainter bursts 
than the PGA or the instruments on XMM-Newton and 
Chandra. Especially for EXO 0748-676 we find a much 
lower SWT fraction than from the XMM and Chandra 
observations. Taking only the PGA bursts into account, 
we find SWT fractions of (11 ± 6)% for EXO 0748-676, 
(13 ± 4)% for 4U 1636-53, and (23 ± 6)% for 2S 1742- 
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294, all of which are within l.Scr from the SWT fractions 
we obtain from Monte Carlo simulations, with an initial 
SWT fraction of 30%. Therefore, in the range of mass 
accretion rates where SWT bursts occur, approximately 
30% of the bursts have a short recurrence time. 

5. CONCLUSIONS 

We studied thermonuclear bursts with short recurrence 
times (SWT) using a large catalog of bursts from multi- 
ple sources observed with the RXTE PCA and the Bep- 
poSAX WFCs. The short recurrence times are of in- 
sufficient duration to accrete the fuel that burns in an 
SWT burst. Bursts are seen to occur in events of up 
to four bursts: double, triple and quadruple bursts. We 
report the shortest recurrence time ever found: 3.8 min- 
utes. We confirm the result of Boirin et al. (2007), that 
SWT bursts are on average less bright, cooler, and less 
energetic than LWT bursts. The decay profiles of SWT 
bursts lack the longer decay component from the rp- 
process, suggesting that SWT bursts take place in a 
hydrogen-depleted layer. 

Some sources exhibit short recurrence times at all val- 
ues of the mass accretion rate where normal bursts oc- 
cur, while for others SWT bursts are limited to a smaller 
range of accretion rate. In this range the fraction of 
bursts with a short recurrence time is consistent with 
30%. Two frequent bursters that likely accrete hydrogen- 



rich matter do not show any SWT bursts. 

Only the hydrogen-accreting neutron stars in our 
catalog exhibit SWT bursts. This suggests that the 
hydrogen-burning processes are responsible for the in- 
complete burning of the available hydrogen and helium 
during bursts. The mechanism for halting the burning 
is still unknown. It will require further theoretical mod- 
eling of hydrogen-accreting neutron star envelopes to re- 
solve this issue. 

As far as we know the spin of the sources with SWT 
bursts, they are all fast rotators. This indicates that ro- 
tational mixing can be responsible for ignition of follow- 
up bursts on a time scale of approxima tely 10 minutes 
(jPiro fc Bildstenl [20071 : iKeek et al.l [20091 ) . The number 
of SWT sources with known spin is small. Measure- 
ments of the spin frequency for more sources and further 
model studies of rotational mixing will better constrain 
this reignition scenario. 
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Fig. 15. — Chandra light curves of EXO 0748-676 as observed in 2001 and 2003. We show the hard and soft energy bands as well as the 
hardness ratio. The observation id is indicated. In the hard band we number the bursts. Figure continues on next page. 
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APPENDIX 

CHANDRA OBSERVATIONS OF EXO 0748-676 

Chandra ACIS-S observations with High-Energy Transmission Grating (HETG) in 2001 and 2003 of EXO 0748-676. 
We use the level 2 event files prepared by the Chandra Data Archive and select events from a circle centered on the 
source and two bands overlapping the signal from the gratings. Using the CIAO software package (version 4.1), we 
extract light curves in the 0.2 — 5keV and 5 — lOkeV energy bands at 60s time resolution, and calculate the hardness 
ratio by taking the ratio of the count rates in the hard and soft bands (Fig. [T5)l . The light curves are similar to the 
XMM-Newton observations performed in 2001 and 2003. Eclipses are present at the binary period of 3.8 h. Especially 
in the data from 2003, dipping is present in the soft band. We, therefore, use the hard band to locate bursts. We find 
a total of 41 bursts; 11 have a short recurrence time. The burst events are divided into 20 singles, 9 doubles, and 1 
triple. The SWT fraction is (27 ± 9)%. Note that burst 22 of observation 4573 appears anomalously long. This is 
caused by a raised detector background level during a few hundred seconds after the burst, combined with the bin size 
of 60 s we used for the figure. 
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Fig. 1151 continued. 



